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Abstract 



The cluster cumulant formula of Kubo is derived by appealing only 
to elementary properties of subsets and binomial coefficients. It is shown 
to be a binomial transform of the grand potential. Extensivity is proven 
without introducing cumulants. A combinatorial inversion is used to re- 
formulate the expansion in the activity to one in occupation probabilities, 
which explicitly control the convergence. The classical virial expansion is 
recovered to third order as an example. 

1 Introduction 

The classical (Ursell-Yvon-Mayer) virial expansion is the traditional introduc- 
tion to real-space particle correlations in basic statistical physics textbooks, 
both old £Q and new PJ. By contrast, cumulant methods, of which the virial 
expansion is a special case, are left to more specialized texts |3J- A possible 
disadvantage of this approach is that physical issues are initially confused with 
purely formal ones, stemming from the specific expression for the grand poten- 
tial il = —kT'fy as a logarithm of a sum: 



The most prominent formal problem is extensivity: it is important to demon- 
strate that the 'factorization property' of matrix elements leads to a f propor- 
tional to the number of single-particle states. 

The formal and physical side of the problem were neatly separated by Kubo 
0]. This article is built around a very compact formal derivation of one of his 
main formulas, appealing only to elementary algebra. Its principal feature is 
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that the usual ordering, by which the cluster cumulant formula is a rearrange- 
ment of the cumulant expansion, is reversed, and the cluster formula appears 
as the basic one. At the level considered here, operative formulas can be devel- 
oped without invoking cumulants at all, physics being determined by the choice 
of expansion parameter. Given the importance of cumulants in general, this 
is not necessarily a pedagogical advantage. It is nevertheless hoped that the 
present work gives a reasonable introduction to the issues involved in cumulant 
expansions, at a level still attainable in a general (graduate) course of statistical 
physics. To this end, the virial series is rederived as an example, not because 
this is the quickest way to do it (it obviously is not), but in order to provide 
a list, in some logical order, of the steps involved in reducing the fully gen- 
eral quantum formula to an operative classical limit. Once such a framework 
is established, each of these steps can be made the starting point of a more 
specialized development. 

Being concerned with interpretation, rather than calculation, this article 
cannot do justice to the numerous, sometimes very refined, applications of the 
cumulant approach, in its broadest sense, which have developed over a pe- 
riod of more than half a century. Outside the original context of a classical 
gas 5 , perhaps the most detailed work was done on spin models, under the 
heading of high-, and low-, temperature expansions 0. Another long-standing 
field of applications are polymers, including fundamental issues, such as the 
excluded- volume problem far from the gelation point . Strong electron corre- 
lations have also been described by cumulant expansions. The common 
denominator of all these applications is the need to treat correlations in real 
space, when neither scaling nor periodicity can be used to simplify the prob- 
lem. An interesting variant of the cumulant method is when a finite system is 
treated in some 'inverse' space, which diagonalizes part of the Hamiltonian; this 
is a classic microscopic approach in nuclear physics |12|. 

Concerning the presentation, I am aware that combinatorial manipulations 
are not part of the usual education of a physicist. If the binomial transform were 
as commonly known as the Fourier transform, the main result in the next section 
would be a 'one-liner'. I have tried to spell everything out, using familiar set- 
theoretic notation. A word of advice to the diligent reader: follow the formulas 
on a small example, with two, or at most three, single-particle states. If it 
suddenly appears easy, that is because it is. 

2 The Kubo formula for fermions 
2.1 The activity expansion 

The diagonal matrix element in TV-particle space is denoted 

U l{N ) = (ii, . . . ,ijv| e~ /3ff \k, . . • ,ijv) , (2) 

where a set of occupied single-particle states (configuration) is denoted by 
i(N) = {ii, . . . , ijv}. It is assumed that the total number of single-particle 
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states is L, and they comprise a set, to be denoted C. In this set-theoretic 
notation, the Massieu potential (JTJ reads 



* = In 



E 



i(N)- 



(3) 



0Ci(AT)C£(L) 



The sum is over all subsets of of the set C of single-particle states, which may of 
course take at most L fermions. It will be assumed throughout that the vacuum 
expectation Uq = 1. All the derivations in this article will refer to fermions. 
Only the final expressions for bosons are given in the Appendix. 

The motivation for the next step is that one would like a sum outside the 
logarithm, not inside it. So define an intermediate expression, 



M 



In 



E 

j(M)C£ $<Zi{N)<Zj(M) 



E 



(AT)- 



(4) 



Here the sum inside is limited to some subset of C, with exactly M states, and 
the sum outside is over all possible choices of such a subset. Notice that \& = *S?l, 
since the sum outside the logarithm then reduces to a single term, j(L) = C. 
Now introduce the pair of inverse relations 



M 



M 

E 

m=0 



rn 



M 



\&(m), 



\£f(m) 



M=0 v 



M- 



(5) 
(6) 



If one of these is inserted into the other, the result is an identity, so they are 
valid independently of what the \& m may be. Taking at M = L, one finds 



(7) 



which is the Kubo cluster cumulant expansion, when ^(m) is given in terms of 
the *m in equation (@}. (Note that *(0) = *o = lnt/0 = 0.) 

To see the structure of Q), look at the first two terms, reverting to standard 
notation for a moment: 




*(1) = 5^1n(H-e^"(i|e- /Jlr |i», 
»=i 

*(2) = Yl In (1 + <i| e"^ H \i) + (j\ e^ m \j)+e 2 ^ 



(8) 



h3\e 



-f3H 



l<i<j<L 



-(L-l)£>(l + e^ (i\e-P H \i)) 
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e^{j\e-^ H \j)+e^{i,j\e-^\U)) 



[ln(l+e^(i|e-^|i) 

\<i<j<L 

- In (1 + (i\ e- fm - In (l + (j| e~ fiH \j))] 

(i,j\e-P H \i,j)-(i\e-P H \i) (j\e-P H \j) 




(1 + (i\ e-? H \i)) (1 + (j\ e~^ H \j)) 



This should be compared with equation 6.22 of reference 0]. The lowest power 
of the activity appearing in ty(m) is e' 3 '™ , so I shall refer to this as the 'activity 
expansion'. If the basis \i) diagonalizes the Hamiltonian, ^(1) gives the exact 
solution of the trivial problem, and the other 'J' (to) are zero. 

Of the three forms in which if? (2) is given, the first is the defining form (JBJ. 
The second is Kubo's original form, where all terms appear under a single sum. 
The original form reads (equation 4.13 of reference @|) 



#(m) = #"("»)]> 



(10) 



i(m)C.£ 



where 



#(m)] = 

0Cj(Af)Ci(m) 



,-M 



hi 



E 

0Cfc(i)Cj(M) 



(11) 



Finally, the third is a compact form, where all logarithms have been collapsed 
to a single one. The trouble with it, although it looks prettiest for small m, is a 
large (hyperexponential) explosion with to in the number of terms in the fraction 
multiplying e^ Mm , as the fraction arises from multiplication of the polynomials 
in which appear under the logarithms in the original form. 

It is easy to prove that the defining and original forms are equivalent. One 
starts from equation (|10f) and simply interchanges the order of summation: 



E E (-i) m - M M-.-j(M)...] 

i(m)C£0Cj(M)Ci(m) 

E (-i) m - M H--i(M)...\ Y. 1 

0Cj(Af)C£ 

Y (-ir- M H...j(M)...} 

0Cj(M)C£ 



j(M)<Zi(m)<ZC 

L - M N 
to - Mj 



= £(-u 



i-Af 



A/>0 



L - M 
to - M 



A/. 



where, in the last step, all M-particle configurations have been grouped together. 
From this point of view, the present article rests on the observation that the 
last line above is formally invertible. 



(9) 
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2.2 Extensivity 

Extensivity itself cannot, of course, be proven without some reference to the 
interactions involved. What will be proven here is more properly called size- 
consistency: if the matrix element E/j( m ) in a configuration i(m) can be expressed 
as a product of lower matrix elements, then the contribution of that configura- 
tion to the grand potential is zero. 

We shall need an elementary property of sets. Namely, if a set i(m) is written 
as the union of two non-overlapping, non-empty subsets, i(m) = ii(mi)Uz2(m2), 
then all the subsets of i(m) may be obtained by writing down all subsets of 
ii(rai) and 12(^2), and combining them in all possible ways. In particular, if 
there is a sum over subsets of i(m), it can be written as two sums: 

]T [...j(n)...}= E [•..ji(ni)Uj 2 (n 2 )...]. 

<l)<Zj(n)<Zi(m) 0Cji(ni)Cii(roi) ®Qh(n 2 )Qi 2 (m 2 ) 

(12) 

This will be referred to below as the 'subset property'. 

The proposition is as follows: let i(m) be the union of two non-overlapping, 
non-empty subsets, i(m) = i\(m\) U 12(^2). Let 

^fc 1 (i 1 )ufc 2 (; 2 ) = U kl ( h )U k2 (i 2 ) (13) 

whenever fci(/i) C ii(m\) and k^li) Q ^O^a). Then tfi[i(rn)] = 0, where 
ip[i(m)] is the contribution of the set i(m) to 9(m) in (|10|l . 

For the proof, first observe that under these assumptions, the sum under the 
logarithm in (|llfl factorizes: 

E u m = E u ^(h) E u w- ( 14 ) 

0Cfc(OCi(n) 0Cfci(ii)Cji(Tn) 9Ck 2 (l 2 )Cj 2 (n 2 ) 

This is trivial: the sum is rewritten by the subset property, and the matrix 
element factorizes by assumption. It follows that the logarithm of l|14|) is the 
sum of two logarithms, one a function of ji(ni) alone, the other of ^2(^2)- The 
contribution of the first logarithm to (|1 1|) reads 

E (-lr-lnl...^)...] (15) 

0Cj(n)Ci(ro) 

E (-l) mi - ni ln[---Ji(»i).--] E (-1)" 12 -" 2 , 

0Cj 1 (m)Ci 1 (mi) 0Cj 2 (n 2 )Ci 2 (m 2 ) 

because the sum here can also be written by the subset property. Since the 
number of subsets of 1(1712) with fixed number n,2 is given by a binomial coeffi- 
cient, 

(-l) ma - B »= ^ f^ ! )(-l) m >- n » = (l-l) m », (16) 

0Cj 2 (n 2 )Ci 2 (m 2 ) n 2 =o\ n2 S 

this is zero for ni2 > 1; similarly the contribution of the second logarithm is 
zero for mi > 1, so the proposition is proved. 
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3 The probability expansion 



In this section, the combinatorial approach is pushed a step further, to rewrite 
Kubo's formula in a particularly transparent way. (The derivations are 
written somewhat more tersely than in the other sections.) Probabilities will 
replace the activity as the expansion parameters, analogously to passing from 
activity to concentration in the classical case. It should be emphasized that this 
only affects the form of the ^[z(m)] in (|1U|I . their value remaining the same, term 
for term, for all i(m) with m > 2. The simplest example of the transformation 
is in the two ways one may write ^(1), 

*(1) = ]T In (1 + e^Ui) = - E In (1 - Pi ) , (17) 

i i 

where the 'occupation probability' 



would be just the Fermi function in the quantum non-interacting case. 



* = YT^k (18) 



3.1 The general transformation 

Here a whole class of ways to expand ^ will be shown to be equal to the activity 
expansion J7J, term for term. In other words, they are merely different ways to 
rearrange the contributions under the logarithms in . A special choice then 
gives the expansion in the probabilities (jT%|) . alluded to above. 

Take an arbitrary set of L variables £j, i = 1, . . . , L. Define the quantities 
W by the pair of inverse relations 

U i(m) = e - " e< J2 W Xn), (19) 

0Cj(n)Ci(m) 

W Kn) = {~l) m - n e + ^^U l[m) . (20) 

$<Zi(m)<Zj(n) 

One could just write Xi instead of e~ l3ei , but the notation is meant to be sug- 
gestive. The grand partition function may now be written 

E e^ N U KN) =f[(l + e^-^) £ (21) 

0Ci(AT)C£(L) i=l <dCj(n)CC(L) 

Here the /i's are just the Fermi functions corresponding to the Obviously, 
they are the new variables of the partition function, replacing the activity. Now 
define, by analogy with equation (0J, 

* M = J2 ln E fh---fj n W j{n) , (22) 

i(M)<ZC $Cj(n)Ci(M) 



G 



and the main statement of this section is as follows: 

$(m) = *(m), m > 2, (23) 

where ^(m) is the binomial transform JBJ of ^m- In other words, the activity 
expansion is unaffected by the transformation, except in the first term; it is easy 
to show that 

E In (l + e^ - *)) + *(1) = *(1), (24) 

i 

with the £, 's cancelling exactly. Of course, since the partition function [left-hand 
side of equation J2U] does not depend on them, they all must cancel in the end; 
but the statement here is that they do so term by term in the expansion, when 
m > 2. 

To prove this, express W in ^>m back in terms of the U's: 
E h-finW m = J] (l + e*"-))" 1 E ^U HN) . 

0Cj(n)Ci(M) iei(M) Q)Ck(N)Ci(M) 

(25) 

Now observe that the sum is the same one as appears in the definition J3J of 
^m, by which ^(m) is defined. The difference ^(m) — ^(m) is thus due to the 
product in front, and reads explicitly, by equation 

Y.^r- M {l~ M M ) E E m(i+e/W), (26) 

M V ' i(M)<ZC iei(M) 

which is, in fact, zero for m > 2. Namely, in the sum over configurations 
(subsets of C), each given single-particle state will appear times, so upon 

exchanging the last two sums, one gets 



L 



i=l 



in (i + £*»-•» 



Y.(-ir-"{tZ M M )( L M -_ 1 ,). <m 



A I 



and after rearranging the binomial coefficients, the last sum is equal to 

i: 1 1 )D-T*C:l)»c-ir', <*) 

' M v 7 

so the statement is proven. 

3.2 The probability expansion 

By choosing 

e-0** = (i\ e~ pH \i) = Ui, (29) 

one finds = pi [equation lfT5)l]. and Wi — 0. This is the useful case, so let us 
denote the W's for this special choice by the letter S, for 'subtracted': 

Si( m ) = E ( — - L ) n_m ^j(n)' (^) 
0Cj(«)C 4 (m) 
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where 



U. 



U 3(n) 



00 = 1. 



In terms of these, the ^[i(m)] in read, for m > 2, 



#(m)] = £ (-ly—ln 

0Cj'(ra)Ci(m) 



Pk 1 ---Pk l S k{i) 

0Cfc(QCj(n) 



noting that 5*0 = 1. This is the 'probability expansion'. For example, 



— 



UiUj 



-1, *(2)= ^ ln(l+ ptpjSy), 



(31) 



(32) 



(33) 



l<i<j<L 

and even \&(3) is short, in the compact form: 



*(3) 



53 In 1 + ptPjPk 

l<i<j<k<L L 

Sijk PiSijSik PjSijSjk PkSikSjk PiPjPkSij SikSjk 



(1 +piPjSij)(l +PiPkS ik ){l +PjPkSjk) 



(34) 



These are higher-order corrections to the non-interacting ('undergraduate') for- 
mula (|17|l . It is obvious how the probabilities control the convergence. 

It should be noted that the 'compact' form has not become really compact, 
but rather that the hyperexponential explosion takes off a little later (because 
Si = 0, so the polynomials being multiplied are shorter). For instance, the nu- 
merator in the compact form of ^(4) has 15 629 terms in the activity expansion, 
and 'only' 505 terms in the probability expansion, still far less practical than 
the 15 distinct terms, 65 additions and 16 logarithms in the original form. 



4 Example 

In this section, the classical virial expansion will be obtained term by term, from 
the formulas developed so far. While such a derivation is nothing new in itself, 
the purpose is to comment on it from the present 'combinatorial' point of view, 
and give a familiar interpretation of the subtracted matrix elements, formally 
introduced in the previous section. 



4.1 The classical limit 

Begin with a Hamiltonian of the form H = K + V, where K is the usual kinetic 
energy, and V a sum of two-body interactions depending on mutual distance. 
The first choice to be made in the formal expansion is, which basis to use for 
the single-particle states. If the momentum basis, which diagonalizes K, is 
used as a starting point, the term Vl/(1) will be the exact solution of the non- 
interacting problem, valid down to zero temperature, and the corrections will 
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correspond to an expansion in quasiparticle occupation probabilities, after a 
canonical transformation to particles and holes. (If this transformation is not 
made, the occupation probabilities of levels below the Fermi energy tend to 
unity, leading to convergence problems.) Such an example, inspired by nuclear 
physics 1 is beyond the scope of the present work. 

If, on the other hand, the position basis is chosen, one is led straight to the 
Mayer expansion. This is particularly easy to see on a lattice. The one-particle 
matrix element is then just the normalized one-particle partition function, in- 
dependently of position: 



Ui 



e -/3(K+V) 



l£exp[-/fe(fc)] = §, (35) 
k 



where L is now the number of lattice sites, and e(fe) is the non-interacting dis- 
persion derived from K. [In the limit of vanishing activity, 11/(1) = Lhi(l + 
e^^Zi/L) — > e^^Zi, which is the classical non-interacting result.] The occupa- 
tion probabilities (|18() are then also independent of position, and to first order 
in the activity, they become equal to the (dimensionless) fugacity: 

eft%/X , M £r ^ (a\ 3 
P 1 + L \X T J ' {6b) 

where the second limit is of high temperature, with a the lattice constant and 
At the thermal wavelength. 

The classical limit for ^(2) is taken in the usual two steps: first, the com- 
mutator [K 7 V] is neglected, being at least of order h, so that disentanglement 
is trivial: 

e -P(K+V) ^ e -/JX e HJV (37) 

after which the two-particle matrix element reads 



le -PK e -0V\ r 



1 e HS«« 



L 2 




Z{~\y cos [k ■ {n - r 3 )] exp [-/fe(fc)] 



.(38) 



Here Vij = V(\ri — rj\), and the interference term from | {n,rj\ki, kj)\ 2 is shown 
explicitly. In physical units k = p/h, so one observes in the second step that as 
% — * 0, it gets 'washed out' by the sum, giving 

*(2)= ln[l+/(e-^«-l)], (39) 

l<i<j<L 

where p is the position- independent occupation probability l|36f) . multiplying 
Mayer's expansion parameter, gij — e~^" iJ — 1. In other words, the subtracted 
matrix elements SVj in the probability expansion become equal to Mayer's pa- 
rameter: Sij — > gij in the classical limit. 
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This is as far as one can go without invoking the probability expansion l|32f> 
explicitly, since only ^(2) is easy to rewrite in probabilities 'by hand'. It is not 
difficult to show that, when Ti — > and all interference terms are neglected, the 
normalized matrix elements <|31[) in the position basis take the familiar form 



where the sum (product) is over all pairs of indices in j(n) = . . . ,j n }- It 
follows that $ is a function of the gij alone, as is well known. The S"s vanish 
with the interaction, as they should, because all the U's in equation H40|) are 
then equal to unity, so they cancel in the definition 1(50) 1 . 

This is a good place to pause, and put the result ((4*0)l in a physical perspec- 
tive. One could have worried: the same one-particle terms (|29l) appear in the 
denominator of the normalized matrix elements l|31[) as in the numerator of the 
probabilities l|18l) . Might they not cancel, leaving the probabilities in expressions 
like (|34() only formally, but not really, in control of the convergence? Not so: 
equation ((4*0)) shows that one should rather expect one-particle ('kinetic') terms 
to cancel between the numerator and denominator of U itself, this cancellation 
being complete in the classical limit. This is just the statement, that in the 
classical limit one can integrate out momenta from the partition function ex- 
actly. Quantum effects are manifested as incomplete cancellation: interactions 
do affect the momentum distribution. However, to believe that a U would not 
tend to a limit when a one-particle matrix element in the denominator became 
small, is to believe that the many-body state in the numerator has an 'infinite 
stopping power' for that particle, if it can avoid making the numerator small 
as well, despite containing that same single-particle state of high momentum. 
Such drastic effects of the interaction on the momentum distribution are not 
unimaginable. In fact, one of the persistent worries in high-temperature super- 
conductivity is that they could preclude any 'semiclassicaP description of the 
conducting electrons [T3J. Nevertheless, it may be (vaguely) concluded, that 
only 'exotic' collective states would spoil the numerical convergence of Kubo's 
expansion for the grand potential, as long as the occupation probabilities l|18fl 
are reasonably small. 

4.2 The virial series 

Going back from probability to fugacity, 



and expanding to second order in /, in the classical limit one obtains 




(40) 




(42) 
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and since it is consistent to write gu = — 1, the term multiplying f 2 may be 
written as the unrestricted double sum 



i=l j=l A 

where it has been used that g^j = g(\ri— Vj\) depends only on the differences, and 
the arrow means the large-volume limit. It is obvious that B\ is just the usual 
'second virial coefficient', except that it is dimensionless; this can be repaired 
by transferring the factor a 3 from the fugacity (j36|) . and putting La 3 = V . 

The third-order result is recovered along the same lines. All that remains 
of ^(3) 134|) to third order in the fugacity is f 3 J2^ijk, where, in the classical 
limit [inserting (gUJ) into (J3DJ], 

Sijk = gij9ikgjk + gij9%k + 9ij9jk + gikgjk, (44) 

from which it is clear that the subtracted matrix elements are the gener- 
ating functions of labelled graphs with n vertices. Adding the terms in f 3 from 
*(1) and *(2), one gets 

'- -'E-'^ • E s «* ] • ( 45 ) 

*<j i<j<k 

To pass from restricted to unrestricted summation, use 

31 E = E =E-E-E-E+ 2 E' («) 

i<j<k i^j^k i,j,k i,j—k i—k 7 j i—j,k i—j—k 

which may be derived by inverting successive expressions of the type 

E = E+E- (^) 

i,3 i^j i=3 

Inserting (|44|l and (|46|) into (|45|l . it becomes 

JE 5 ^/ 3 ' ( 48 ) 

so that, just as in second order, all that the lower terms in (|45|l do is to remove 
restrictions on the sum in the highest one. In the large-volume limit, the first 
term in (|44|l gives 

\ Es«0ifc$j* ~> \ E 3(Ai)g(A 2 ) 5 (|A 1 - A 2 |) = LB 2 , (49) 

i,3,k Ai,A 2 

and the remaining three give a total contribution 

\ E 9ij9ik - ~ E E 3( A i)5(A 2 ) = 2LBf, (50) 

i,j,k Ai A 2 
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because the g's depend only on the differences, so one finds, finally, 



* = Lf + LB.f + L(2Bf + B 2 )f + 0(/ 4 ), 



(51) 



which is the virial scries to third order, in one notation 

The statement that the interacting problem is 'reduced to quadrature' in 
the classical limit is interpreted here, that the coefficient of f n in \I/ becomes an 
unrestricted sum, namely 



From this point of view, in the classical limit there appears a 'conspiracy of 
terms' which removes quantum restrictions from the sums in Kubo's formula. 

5 Discussion 

This article gives a combinatorial interpretation of Kubo's cluster cumulant 
expansion, as a binomial transform of the grand potential. It primarily explores 
the pedagogical implications of having such a short, but formal, derivation. The 
idea is to develop a self-contained, general point of departure to treat problems 
which require a formulation in real space, assuming only undergraduate prior 
knowledge. In particular, it is found there is no need to introduce cumulants 
explicitly at this level, in order to produce operative size-consistent expressions. 
This was demonstrated in detail for the classical limit. 

In this approach, the distinction is kept between rearranging the series for 
which (effectively) resums different infinite classes of terms, and generating the 
terms themselves, which requires the evaluation of matrix elements. It becomes 
clear in principle, how different choices of expansion parameter (probability, fu- 
gacity, coupling constant,. . . ) necessarily yield different rules for which terms 
appear, and how either quantum entanglement, or topological restrictions from 
the Hamiltonian, both of which change the form of (|4U[) . can spoil the 'conspir- 
acy' by which the classical expressions simplify. 

The basic operational problem in quantum mechanics is to replace sums by 
integrals, or, in more sophisticated language, to pass from functions defined on 
sets (of quantum states) to functions of real numbers (parameters of the Hamil- 
tonian). This is trivial with unrestricted sums, which is the combinatorial con- 
tent of the simplification in the classical limit. The converse is formally the most 
difficult problem of strong correlations: when dynamical effects restrict a (mul- 
tiple) sum to an 'arbitrary' subset of discrete states, i.e. such that no ordering 
can be defined on it, there is no controlled way to express the sum as an integral 
in the large volume limit. Otherwise, the standard way to obtain integrals is 
to introduce ordering by the time variable, leading to Feynman diagrams. The 
relationship between the cluster cumulant expansion and the diagrammatic ap- 
proach has been discussed by Dunn , for the case of a particular self-energy. 
He showed that cutting the expansion off at m-th order amounted to calculating 




{coefficient of f n in . 



(52) 
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all diagrams with at most m-fold momentum integrals exactly, and all others 
approximately [the logarithms in \l/(m) necessarily generate diagrams to infinite 
order] . 

On the other hand, rewriting Kubo's formula as a probability expansion, 
equation (|32|) . shows that convergence can be expected even when one does 
not have the complete solution of the problem. All that is required is that the 
probabilities in the chosen basis are bounded away from unity (and zero). The 
prototype for this is precisely the real-space basis, because position states are 
never stationary, due to the uncertainty principle. 

Two other properties of the cluster expansion are readily obtained. First, it 
was stressed by Kubo @] that all described operations remain exactly correct 
even if the 'matrix elements' are not c-numbers, and the exponential functions 
are replaced by various rules. This follows directly from the fact that the bino- 
mial inversion is a formal identity. For example, when the probability expansion 
can be written in Fock-space operators, the Pauli principle is exactly preserved, 
even if one stops at first order [§]. Second, if matrix elements are used, the 
Pauli principle cannot be satisfied for the whole assembly of particles, as soon 
as the expansion is cut off. Its form then indicates that antisymmetrization is 
taken into account by an 'inclusion-exclusion' procedure; for instance, the same 
two-particle terms, involving Sy, appear with opposite sign in '5(2) and ^(3) 
[equation So the cluster cumulant series is expected to alternate, whenever 

many-body correlations are important. This can be striking in practice |15|. 

To conclude, the second and fourth sections of this article give a compact 
and hopefully useful introduction to an established general treatment of cor- 
relations in real space. In the third section, it is shown that the parameters 
controlling convergence can be interpreted as probabilities, and their associ- 
ated subtracted matrix elements appear as the basic building blocks of more 
elaborate calculations. 
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A The expansion for bosons 

The probability expansion l|32|) is formally the same for fermions and bosons, 
only the definitions of the various quantities change. For bosons, the occupation 
probabilities are 

where i k means, i-th. state, occupied by k bosons. In the non-interacting case, 
this reduces to the familiar e' 3 ^ _£i \ justifying the use of the term 'occupation 
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probability'. 

The normalized matrix elements (|31ll become, for bosons, 



V°° P (k 1 +...+k n )Pn /xfci ,-fc„ 

Z J fc 1 ,...,fc„=l e \J1 •••Jn 


e -/3H 


,-fcl 

Jl • 










./"))•••(>:: : 




(in™ 







(54) 

To lowest order in the activity, these expressions are of course equal to the 
fermion ones, which accounts for the classical limit. Note, finally, that the 
denominator of Uji n ) contains the 'occupation numbers' 

oo 

rn = 5^ e*"" <t*| e-^ |**> , (55) 
fe=i 

in terms of which pi — ni/(rii + 1), as with non- interacting bosons. 
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